Non equilibrium statistical physics with a fictitious time 
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Problems in non equilibrium statistical physics are characterized by the absence of a fluctuation 
dissipation theorem. The usual analytic route for treating these vast class of problems is to use 
response fields in addition to the real fields that are pertinent to a given problem. This line of 
argument was introduced by Martin, Siggia and Rose. We show that instead of using the response 
field, one can, following the stochastic quantization of Parisi and Wu, introduce a fictitious time. In 
this extra dimension a fluctuation dissipation theorem is built in and provides a different outlook 
to problems in non equilibrium statistical physics. 

PACS numbers: 05.10.Ge 
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The fluctuation dissipation theorem (FDT) is the cen- 
tral feature in the study of dynamics of statistical sys- 
tems, when the deviation from equilibrium is small. The 
response of such a system to a small external probe is 
studied within a linear response theory. This is a macro- 
scopic effect which can be observed in a dynamic mea- 
surement. The FDT relates this response to the fluctua- 
tion properties of the thermodynamic system. In its sym- 
plest form it appears in Brownian motion-the molecular 
fluctuation responsible for the agitation of the suspended 
particles are the same as those responsible for the macro- 
scopic friction experienced by the particles in motion. 
Generalised, it says that in an equilibrated system under 
a small perturbation there is a balance between a system- 
atic force and a chaotic fluctuating force. The theorem is 
valid for all dynamical systems in which an equilibrium 
distribution is reached at infinitely long times. Clearly, 
for a system which is far from equilibrium this theorem 
is absent. In this communication, we would like to point 
out that following the technique of stochastic quantiza- 
tion, we can introduce an extra time like dimension in 
these far from equilibrium problems and in this fictitious 
time introduce a FDT. This greatly simplifies the study 
of possible scaling solutions in the problem. 

The FDT facilitates the study of dynamics near equi- 
librium by the existence of the fluctuation dissipation 
theorem which relates the correlation function to the re- 
sponse function. One of the primary difficulties in the 
study of dynamics far from equilibrium is the absence 
of a fluctuation dissipation theorem leading to indepen- 
dent diagrammatic expansions for the response function 
and the correlation function. The difficulty stems from 
the fact that the equilibrium distribution is not known 
and the only way in which we can do an averaging is 
the averaging over noise. The spectacular success in the 



last few decades in the study of dynamics has been in 
the near equilibrium dynamics near second order phase 
transitions. Success in this case implies a close and de- 
tailed correspondence between theory and experiment in 
a variety of systems 0, 0> EL IS E3- The dynamics 
in this case has typically been that of a n-component 
field 4> a (ot = 1, 2, n) satisfying an equation of motion 
(model A, model B etc.) 



Of 



= -r- 



SF 



(1) 



where F is a free energy functional (most often the 
Ginzburg-Landau free energy), with the Gaussian white 
noise rj(f, t) satisfying 



t 2 ) 
(2) 

It is straightforward to check by going to the associ- 
ated Fokker Planck equation that exp(—F/kT) is indeed 
the equilibrium distribution. Since the dynamics is al- 
ways close to equilibrium, it is the averaging with the 
distribution exp(—F/kT) which one has in mind. There 
have been more complicated equations of motion with 
the structure 



— = Va {[4> p ])-T— 



(3) 
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(e.g. models E, F, G, H, J of dynamic critical phenom- 
ena), but the equilibrium distribution has been main- 
tained because ^f- = V a ([<j)p]) keeps the free energy F 
constant in time (in fact that is the motivation behind 
the construction of V a ). 

Models of non equilibrium statistical mechanics do 
not evolve towards an equilibrium distribution. Conse- 
quently, the averaging involved in the construction of cor- 
relation function have to be done explicitly over the noise 
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inherent in the statistical dynamics (the noise written 
explicitly in Eqs QJ and (J5J). One of the most-studied 
models in this category is the Kardar Parisi Zhang (KPZ) 
equation, where the dynamics ( physically the growth of 
a surface by deposition of atoms on a substrate ) is given 
by 



dcj>{r,t) 
dt 



vV 2 (j)~ ^(W) 2 +7] (4) 



(0) 



Integrating over rj, Z = f Ufiexp — S((j)), with the ac- 
tion 



where 



< r ) (n,t 1 )7i(fi,t 2 ) >= 2D S(n - r 2 )5{h - t 2 ) (5) 



The KPZ equation models a growing interface which 
develops from a random deposition of atoms on a sub- 
strate of dimensions D. The randomness in the deposi- 
tion process brings in the noise term in the time evolution 
of the height of the interface. The deterministic part of 
the evolution coming from the deposition and evapora- 
tion of atoms is governed by the chemical potential /i 
which is determined by the total height gradient and its 
derivatives. Clearly a linear term in the gradient is not 
acceptable since the evolution cannot depend of the sign 
of the gradient. Thus the linear term in \i has to be the 
second derivative of the height and the non linear term is 
the square of the gradient. Without the non linear term, 
one has the Edwards Wilkinson model 8] of growth and 
with it the KPZ model. One of the primary uses in the 
models of growth is the question about the roughness of 
the interface. This is expressed in terms of the correlation 
function < <fi(x + r,t)<fi(x,t) >, which scales as r 2a for 
large r. If a > 0,then the height fluctuations are coupled 
at long distances at the surface is rough. If a < 0, the 
fluctuations decay and the surface is smooth. The EW 
model has a = (2 — D)/2 leading to a rough interface 
for D < 2. For the KPZ model, a > and the surface 
is always rough for D < D c , an upper critical dimen- 
sion, while for D > D c , the surface is smooth for small 
values of the coupling constant . The value of D c is still 
controversial @, M,M, El 

Now Eq. igj is not of the form of Eq.QJ (vV 2 (j) ~ 
■|(V0) 2 cannot be written in the form SF/S(f>). It is of 
the form of Eq. ©, with F = {v/2) J d D x(V(j)) 2 , but 
identifying V^([0]) as — (A/2)(V0) 2 , it can be shown that 
F is preserved only in D = 1 (D is the substrate dimen- 
sion). 

The field theoretic technique of handling such problems 
was first explained by Martin, Siggia and Rose an d 
perfected over the years by Janssen Bausch et al 
De Dominicis and Peliti Frey and Tauber [l7j . 
Plischke et al 01 an d a host of other researchers. One 
begins by writing the "partition function" 



5(0) 



D J [ dt 



(J) 



as given by Zee. 

It is with the action that all correlations of the form 
< 0(ri, ti)<f)(f2, fa) > have to be determined. The re- 
sponse functions are introduced by writing Z in terms of 
an auxilliary field, so that 



(8) 



Z = J V[(/)]V[4>]e- s ^ 
where the new action is given by 



Ste,$) = J dtj d D f[0(r,t)[0-^V 2 0--(V0) 2 ]-^ 2 ] 

The way of handling this case of non equilibrium sta- 
tistical mechanics delineated clearly in Frey and Tauber 
p~7| . Our proposal, here, is to use the action of and 
exploit the principle of stochastic quantization to set up 
an alternative approach to the study of this class of prob- 
lems. Stochastic quantization is a method of quantisation 
proposed by Parisi and Wu 0, |2£j based on stochastic 
Langevin dynamics of a physical system in a fifth time 
r. They showed that at the perturbative level, the usual 
quantum field theory was recovered in the limit r — > oo 
of this dynamics. Euclidean quantum field theory corre- 
lation functions for a field 4> corresponding to an action 
S{4>) are given by 



< | T^{x x )4>{x 2 ) <j>{xi) | > 



(10) 



Parisi and Wu proposed the following alternative 
method: (a) Introduce an extra fictitious "time" r in 
addition to the four space time and postulate a 
Langevin dynamics 



d(j)(x, t, t) 
dr 



SS 

"=T+f(.X,t,T) 



(11) 
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where / is a Gaussian random variable with 
< f{x,t,T)f(x',t'T') >= 2S(x-x')S(t-t')S(T-r') (12) 

(b) Evaluate the stochastic average of the 
fields cj>{x,t,T) satisfying Ea.Qll|) ie evaluate 

< 4>(xi,Tl)(f>(x2,T2)...(t>(xi,Tl) >n 

(c) Set T\ = T2 = ■■■■= ti = t and take the limit 
t — > oo. Then one has 

lim < (/)(xi,t)4>(x2,t)....<I)(xi,t) > v 
J P[0]0(x 1 )0(x 2 )....0( a ; ; )e~ s ^ 

(13) 

We now return to the KPZ equation and write the 
action of Eq. in momentum-frequency space as 

s = m I + »'*<-'"• -> 

+ — y^(— iw + vk 2 )q.(k + q)(f>(k, w)(j)(q, w') 

<j){—k — q, —w — «/) — — + vk 2 )q.(k — p) 

4>(q, w')4>{—k — w)4>{k — q,w — w ) 
"T £ [p-{k-p)U(k + q)] 

p,q,w' ,w" 

(j)(p, w')(j)(k — p,w — w f )(j){—k — q. —w — w ,, )4>(q 7 w f/ )} 

(14) 

in accordance with standard results [2l|. The Langevin 
equation in the fifth time r is now written as 



dr 



SS 



6<P(-l,-n,r) 



+ f 



(15) 



with 



< f(i, n, r)/(F, n', t') >= 2S(T+ V, n+ ft', t - /) (16) 

After straightforward algebra Ea. Q15f) acquires the 
form 



8t 



2D Q 



-^^[2(-w + rf)l?- fc) 



4D, 



o r 



- (in + vl 2 )k.(l - k)]cj){k, w, r)(f>(l ~k,n-w,r) 
+ 45- E [^(r-fe)][g.(fc-9)]^-;^,r) 

<£(£ -q,w- w',T')(j){T- k,Q-w, t') + /(i, n, t) 

(17) 



For the free field theory, (Edwards Wilkinson model) 
A = and we get the response function of the system as 



G = (-iw T 



w 2 + v 2 k A ! 
2D 



(18) 



where w T is the Fourier transform variable corresponding 
to the fictitious time r. The correlation function is 



C (k,w,w T ) = [w 2 + 



( W 2 + V 2 k i ) 2 j 



4Z? 2 



-7mG (19) 



as required by the fluctuation dissipation theorem in 
the fictitious time. To get the correlation function of 
the original theory, we need to consider C{k, W, Ti,T2); 
set T\ = T2 = t and let r — * oo. The part which 
survives when t — > oo, is obtained directly from the 
equal-r from Eq. l|19fl . ie. the result which is obtained 
by integrating the right hand side of Ea. l|19l) over w T . 
This yields 2Dq(w 2 + j/ 2 /c 4 ) -1 for the correlation func- 
tion of the Edwards Wilkinson model, leading to the re- 
sult a= (2 - D)/2. 

With A / 0, we develop as usual the fully dressed 
Greens function from perturbation theory and write 



G 



w 



2 1,4 



V Z k 



~IW T 



Dn 



(20) 



The point of our approach is that the w = part of £ 
gives the flow v and the k — part of £ gives the flow 
of D. Further, mode coupling analysis will directly yield 
the exponent z from a power counting analysis. As in 
other studies the result a + z = 2 holds in this case as 
well. 

If we inserted in £(fc, w, w T ) to 0(A 2 ), then we note 
that the contribution comes from two different sources - a 
one loop contribution from the 0(A 2 ) term in Ea. (|17fl and 
a two loop contribution from the O(A) term. The contri- 
bution of the 0(A 2 ) term can be read off from Eq. ((T7|l by 
contracting two of the ^-fields. The resulting correction 
to (n 2 + v 2 I a )/2Dq does not have any new momentum 
dependence and hence it is the second order contribu- 
tion from the 0(A) term which is of significance. The 
expression for £(?, n, w T ) is found after standard alegbra 
as 



A 2 



dPk dw dw T 



e(l n, w T ) = ., , 
V ' ' ' 8D 2 J (2ir) D 2vr 2vr 

{[2(-iw + vk 2 )l- (I- k) - (-n + vl 2 )k ■ (I- k)] 

[2(-iw + in + u(l- k) 2 l- k) + (iw + vk 2 )l- (I- k)] 

G(k, w, w' T )C(l — k,Q — w,w T — w' T ) 

+similar term interchange of k and (I — k) 

except in the first factor} 

(21) 



4 



We now note that the Greens' function can be written 
to 0(A 2 ) as 

G~\l, w, w T ) = -iw T + -^-(n 2 + iA y Z 4 ) (22) 
where 

v% ff l* = v 2 l A + Y;(l,n,w T ) (23) 

or 

u efs l 2 ~vl 2 + ^Z(l,fL,w T ) (24) 

or 

Aul 2 = -±-Z(l,n,w T ) (25) 

In a self consistent mode coupling, we now replace v by 
Av in Eq.J2U, use G as given by Eq. l|22l) and C as follows 
from the fluctuation dissipation theorem. We can carry 
out the momentum count of Eg. I|21|l .keeping in mind that 
vl 2 ~ l z , to find 

£ ~ l D l z l 2z+4 r 4z = l D+4 - z (26) 



Using this in Eq. (|2*S|l . we have l z ~ l D+2 z , leading to 



z = 1 + - (27) 



Our analysis is not valid for D > 2, because then £ 
no longer dominates I 2 . We have checked to that part 
from Eq.J2U and Ea. (|26|l . one can write down flows for 
v Dq and as in the usual analysis, there is a critical D, 
1 < D < 2, where the flow of the coupling constant 
X 2 Dq/i' 3 changes sign. We do not consider that to be 
of any particular significance one way or another. Our 
primary claim is that by introducing a fictitious time, we 
are able to bring back to fluctuation dissipation theorem 
and this allows a power counting scaling analysis for the 
KPZ system which had not seen possible before. 

In closing, we would like to reiterate that the standard 
practice in dealing with out of equilibrium situations has 
been to use a set of fictitious fields called response fields, 
which provides a field theoretic prescription for the re- 
sponse function. Instead, we have proposed the intro- 
duction of a fictitious time in which a FDT holds and 
thereby only correlation functions need to be calculated. 
Study of scaling solutions is thereby greatly simplified. 
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